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We consider the possibility of topological quantum phase transitions of ultracold fermions in 
optical lattices, which can be studied as a function of interaction strength or atomic filling fac- 
tor (density). The phase transitions are connected to the topology of the elementary excitation 
spectrum, and occur only for non-zero angular momentum pairing (p-wave, d-wave and f-wave), 
while they are absent for s-wave. We construct phase diagrams for the specific example of highly 
anisotropic optical lattices, where the proposed topological phase transitions are most pronounced. 
To characterize the existence of these topological transitions, we calculate several measurable quan- 
tities including momentum distribution, quasi-particle excitation spectrum, atomic compressibility, 
superfluid density, and sound velocities. 



PACS numbers: 05.30.Fk, 03.75.Hh, 03.75.Ss 

In the last few years there has been a tremendous in- 
terest in atoms loaded into optical lattices, where many 
known phases of standard condensed matter physics can 
be simulated under a tighly controlled environment. Ex- 
amples of such realizations are the observation of the 
superfluid-to-insulator [H transition in atomic Bose sys- 
tems, metallic molecular [§] and superfluid [H phases 
of atomic Fermi systems. The production of quantum 
degeneracy of atoms and molecules in optical lattices has 
allowed the merger of standard condensed matter and 
atomic and molecular physics into a vibrant field of re- 
search called condensed atomic and molecular physics. 

In standard fermionic condensed matter systems it has 
been very difficult to study systematically the effects of 
strong correlations as a function of particle density and 
the ratio of hopping to interaction strength, as the abil- 
ity to control these parameters is generally very limited. 
However, in atomic systems, the strength of atom-atom 
interactions, the atom hopping, and the atomic filling fac- 
tor in optical lattices are easily controllable and allow the 
exploration of a much wider phase space. In particular, 
the ability to control interactions relies on the existence 
of Feshbach resonances, which have been used to study 
fermions loaded into optical lattices both in the s-wave 
channel 0] for 6 Li, and the p-wave channel [f| for 40 K. 
More recently, p-wave Feshbach resonances in 6 Li were 
used to produce p-wave molecules [1, 0] in harmonically 
trapped systems, and it is likely that similar experiments 
may be attempted in optical lattices. 

Currently, the only fully confirmed fermion superfluid 
with pairing higher than s-wave, is liquid 3 He, which is 
not found on a lattice. Very few lattice condensed mat- 
ter candidates exist, such as the ruthenates [8(, and or- 
ganics [9(, but since the control of carrier density, hop- 
ping and interactions in these systems is very limited, the 
exploration of the phase diagram of alleged p-wave lat- 



tice superfluids has been hindered. Thus, the search of 
higher-angular momentum superfluidity in atomic sys- 
tems is very important not only to help ellucidate the 
symmetry of the order parameter in analogous condensed 
matter systems, but also to explore new phases that are 
not accessible in standard condensed matter due to the 
lack of control of interactions, density and hopping. 

Since Feshbach resonances in optical lattices have al- 
ready been observed in the p-wave channel @, the in- 
teraction strength can be tuned continuously from weak 
(BCS) to strong (BEC) attraction limits. Here, we 
take advantage of the tunability in optical lattices to 
study theoretically the occurance of unusual topologi- 
cal quantum phase transitions in three-dimensional but 
anisotropic optical lattices with superfluid order parame- 
ters in non-s-wave channels. We construct the phase dia- 
gram for anisotropic three-dimensional optical lattices in 
the fermion density versus interaction strength plane and 
identify up to five different quantum phases for non-zero 
angular momentum states depending on the momentum 
space topology of the quasiparticle excitation spectrum. 
For some p-wave states, we find that the quasiparticle 
excitations are gapless in the BCS regime, and are fully 
gapped in the BEC regime. To characterize the change 
in topology of the quasiparticle excitations, we show that 
the momentum distribution, atomic compressibility, su- 
perfluid density and sound velocity are non-analytic func- 
tions of the interaction strength exactly where the topo- 
logical changes occur. 

Hamiltonian: To describe the physics described 
above, we study ultracold fermions in anisotropic three- 
dimensional optical lattices described by the single band 
dispersion 

£k = — 2t x cos(k x a) — 2t y cos(k y a) — 2t z cos(fc z a), (1) 

Here, the hoppings t x > t y > t z are chosen to be differ- 
ent and a is the optical lattice spacing. We work with 
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the Hamiltonian H = Hkin + Hint , where the kinetic en- 
ergy part is H kin = J2k,a ^t^K^ with £k = «k - Mi 
where the chemical potential /i may contain the standard 
Hartree shift. The interaction part of the Hamiltonian is 

H int -|EE Vdhsfr k ')^(k, q)6 75 (k', q) (2) 

kk'q 0.(3^8 

with 6^(k, q) = V'l k+q/2ia V'k+q/2,^ where the labels 
/?, 7 and S are the pseudo-spin indices and the labels k, 
k' and q represent linear momenta. We use units where 
h = ks = 1, and allow the pseudo-spins indices to take 
two values (pseudo-spin S = 1/2) corresponding to two 
hyperfme states labeled as | f) and | |). 

In the case, where the hyperfme states (pseudo-spin 
indices) and the center of mass coordinates are uncoupled 
the model interaction tensor can be chosen to be 

Va^s^M) = -Vr<Mk)^r(k')r Q/ 3 7 ,5, (3) 

where the tensor T a f3 1 s — 8 a p8^/2, for the pseudo- 
singlet pairing case (5i,2 = S\ + S% = 0); and where 
r Q/ 3 71 5 = v Q/3 -v 7i5 /2 with v a p = (i<T<jy) a p, for the pseudo- 
triplet pairing case (Si ,2 = Si + S 2 = 1). Here, V? has 
dimensions of energy and represents a given symmetry of 
the order parameter with basis function 0r(k) and r (k') 
representative of the orthorhombic group (Z?2/i)- 

Self-Consistent Equations: At the saddle point, the 
pairing field V\(ki -fk 2 , r) is taken to be r independent, 
and to have center of mass momentum ki + k 2 = 0, 
becoming £> (ki+k 2 ,r) = A r <5 kl+k2 , + 5V (k 1 +k 2 , r) 
for the singlet case, and 2?i(ki +k 2 , r) = ^Ar^k 1 +k 2 ,o + 
<52?i(k! + k 2 ,T) for the triplet case. The corresponding 
order parameter equation is 

1 = J2 l^r||0r(k)| 2 tanh(/?£ k /2)/2£; k , (4) 

k 

The number equation is obtained from N = —dQ/d[i, 
where (3fl = — In Z is the thermodynamic potential and 
Z = Tr(cxp —f3H) is the partition function, leading to 

N = N Q + iVflact, (5) 

where N = X) k ™k, and n k = [1 - £ k tanh(/3£ k /2)/£' k ] 
is the momentum distribution. The additional term 
-'Vfluct = — 3f2fl uc t/9/x, where Ofl uct are Gaussian fluctua- 
tions to saddle point f2 - These two equations must be 
solved self-consistently in order to provide the order pa- 
rameter amplitude Ar, the chemical potential fi, and the 
quasiparticle excitation energy 

E* = yJ& + \A r \*\MV)\ 2 - 

Order Parameter Symmetries: For singlet pairing the 
saddle point field is = Ar, while for the triplet case 
it is = r/i Ar, which is related to the d- vector order 



parameter by di(k) = Ylk^i Vr(k). For orthorhombic 
lattices without breaking time-reversal or parity and no 
coupling between pseudospins and center of mass coordi- 
nates, the only order parameters for superfluidity allowed 
by symmetry are: (a) A(k) = Ar0r(k) for singlet states 
and (b) d(k) = f/Ar</>r(k) for triplet states. This im- 
plies that there are only 8 symmetries allowed for the or- 
der parameter which are consistent with the orthorhom- 
bic D2h group. There are four options for the singlet 
case: (a) s-wave with A(k) = A s ; (b) ci^-wave with 
A(k) = A dwV XY; (c) d xz -wave with A(k) = A dxZ XZ; 
and (b) d yz -wave with A(k) = A dyZ YZ. There are also 
four options for the triplet case: the d-vector in momen- 
tum space for unitary triplet states in the weak spin- 
orbit coupling limit is characterized by one of the four 
states: (a) p^-wave with d(k) = r)A Px X] (b) p^-wave 
with d(k) = 7]Ap y Y; (c) p z -wave with d(k) = t)A Pz Z\ 
and (d) f xyz -wave with d(k) = f]A^ y ^XY Z . Since, 
the Fermi surface can touch the Brillouin zone bound- 
aries the functions X, Y, and Z need to be periodic and 
can be chosen to be X — sin^a), Y = sin(k y a), and 
Z = sin (k z a). The unit vector r) defines the direction of 
d(k). From here on we scale all energies by i^and choose 
t y /t x = 0.2 and t z /t x = 0.008 such that t x > t y > t z 
where the largest number of non-trivial phases occur. 

Topological Transitions: We discuss three distinct 
phases based on the normal state and quasiparticle Fermi 
surface (FS) topologies and the nodal structure of the or- 
der parameter as a function of /j,. We consider the "nor- 
mal state" FS defined in the first Brillouin zone by = 0, 
keeping in mind the periodicity in fc-space, and define the 
special values /i* = —2t x — 2t y — 2t z : = —2t x — 2t y +2t z ; 
^3 = -2t x + 2t y - 2t z ; n% = -2t x + 2t y + 2t z for filling 
factors < N < 1. The order parameter always has no 
nodes for the s symmetry, while the nodes for the pi sym- 
metry can only occur on the planes ki = 0, ±7r where i 
is x, y, or z. For the superconducting state, the intersec- 
tion of the FS and order parameter nodes constitute the 
loci of gapless excitations, where E(k) — 0. Since E(k) 
is always gapped for the s-wave symmetry, there are no 
quantum phase transitions present. However, topological 
quantum phase transitions occur for pi symmetries. 

For all triplet p^-wave symmetries, E(k) is fully gapped 
only for /i < since there is no Fermi surface. For 
fil < /i < the order parameter nodes intersect the 
Fermi surface and hence quasiparticle excitations are 
gapless. For fj,± < // the Fermi surface splits into two 
sheets that separate along k x so that the p x nodes no 
longer intersect the Fermi surface opening a gap in the 
excitation spectrum. However, the p y , and p z nodes still 
intersect the Fermi surface so excitations remain gapless. 

Representative phase diagrams for the p y symmetry 
are shown in Fig[T] There are three distinct superfluid 
phases characterized by Fermi surface connectivity and 
quasiparticle excitation spectrum: (1) no Fermi surface 
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FIG. 1: a) The T = phase diagram of filling factor N ver- 
sus interaction V Py . The cross indicates the location of a 
suspected p-wave superfluid in highly anisotropic lattices of 
standard condensed matter, see Ref. [IJ. b) the tempera- 
ture versus interaction phase diagram at N — 0.5 (quarter 
filling) for the p y symmetry. The solid lines indicate asymp- 
totic forms of the critical temperatures for weak and strong 
coupling, and the dotted lines indicate the phase boundaries 
between topological phases of the type (1), (2) and (3). The 
hopping ratios used are t y /t x = 0.2, and t z /t y = 0.008. 



and fully gapped for all pi symmetries (fi < //*), (2) 
one sheet Fermi surface and gapless for all pi symmetries 
(pL\ < fx < nX)-> (3a) two sheet Fermi surface and gapless 
for p y , p z (pL% < fi) (3b) two sheet Fermi surface and fully 
gapped for p x (/I4 < /1). In addition, phase (2) splits into 
three regions using the finer classification of Fermi surface 
topological genus: (2i) genus zero (fi\ < [i < pL^), (2ii) 
genus one < < /X3), (2iii) genus two (fit, < \x < /U3). 
There are several qualitative features of interest in the 
phase diagrams. For a fixed low density and all triplet 
symmetries, the chemical potential does not go below the 
bottom of the band until a critical coupling is reached 
indicating the formation of a bound state requires a finite 
interaction strength, which occurs at \V P /2t x \ = 3.5052, 
\V f 

Py 



/2t x \ = 2.3952, and \V p j2t x \ = 1.8150, for the p x , 
, and p z symmetries respectively. 
Momentum Distribution: While changes for s-wave are 
smooth, dramatic rearrangements of the momentum dis- 




FIG. 2: Momentum distribution in the plane k z = for the p y 
symmetry for = 0.5 and in the vicinity of \V Py /2t x \ ~ 11 
when (a) < /x* and (b) > fj,*. Note the jump in the 
momentum distribution at the origin from a) to b), signaling 
a topological transition at /i = /i^. 



tribution n(k) accompany the topological transitions for 
p-wave d-wave and /-wave. For brevity, we analyze only 
the pi symmetries where i — x,y, z. Near the points 
where the Fermi surface changes topology (connectivity 
or genus) n(k) » 1 + sgn(fy) -sgn(r5 M ) A 2 pj 1^- 1 2 /(2|^| 2 ) 



and <5k = k 



k* where j 



1,2,3,4. 



where 5/j, = /i — /1* 
The momentum k* is the point (up to lattice periodicity) 
where the Fermi surface changes topology: k* = (0, 0, 0), 
k^ = (0,0,tt), k^ = (0,tt,0), k| = (0,7r,7r). This ex- 
pansion clearly breaks down for dfi — 0, where we find 
n(k) w 1 - \tj\\5ki\/Ao for |«5k,-| 0, and n(k) w for 
|r5kj| = 0. When /i crosses the boundaries there are 
clear discontinuities in n(k). Plots of n(k) for the p y 
symmetry are shown in Fig. [2] for N = 0.5 and in the 
vicinity of ^ w p* (\V Pv /2t x \ « 11). 

Compressibility and Superfluid density: The atomic 
compressibility and superfluid density provide additional 
signatures of topological quantum phase transitions in 
optical lattices. These quantities can be directly ex- 
tracted from an analysis of quantum (T = 0) phase fluc- 



tuations. The effective action S e ff 



7(0) 



SS ' c ff has a 



saddle point S^J = /30o and fluctuation rS^cft = (3flauct 
contributions. Integrating the amplitude fluctuations in 
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FIG. 3: Plots of (a) atomic compressibility k and (b) sound 
velocity c y versus interaction V Vy for the p^-symmetry at zero 
temperature. The parameters are the same as in Fig. \T\ 



SS c ff leads to the phase only action 



SS, 



phase 



= I E t A K) 2 + PijW) 0{q)e{-q), (6) 



with 6(q) being the phase field of the order parameter and 
A = N 2 K /V, where « = -i (f£) riV = ^ (f)^ is 

the atomic compressibility and pij — V^ 1 ^ k [n^didj^] 
is the superfluid density tensor. Notice that AS does 
not only give k and p^ but also the phase-only collec- 
tive mode (sound) velocity w(k) = 



C X Q X 



-yVy 



where = p xx /A, c 2 y = p yy /A and c\ = when 
the analytic continuation icj„ — > cj + id is performed. 

Measurement of two quantities out of three (sound ve- 
locity, superfluid density and compressibility) can yield 
the third one. In particular, the techniques used to mea- 
sure accurately the sound velocity of superfluid fermions 
in harmonic traps 0, [n| and to measure compressibil- 
ity [IH can be used in the lattice case. In Fig. [3l we show 
(a) the atomic compressibility k, and (b) the sound veloc- 
ity c y — y/ (p yy V)/(N 2 K,) along the y-direction. Notice 
the clear non-analytic behavior in these two quantities 
when the phase boundaries between different topological 
phases are crossed. 

The non-analytic behavior of the momentum distribu- 
tion, sound velocity and compressibility, as function of 
interaction strength or density (not shown) character- 



izes the topological quantum phase transition described. 
Since topological transitions occur without a change in 
symmetry, they can not be classified under Landau's 
scheme, however discontinuities occur in the third deriva- 
tive of the thermodynamic potential (first derivative of 
compressibility), and the transition can be identified as 
a third-order transition under Ehrenfest's classification. 
This type of topological quantum phase transition be- 
longs to a class of transitions first envisioned by Lif- 
shitz [l3| for non-interacting Fermi systems, where the 
Fermi surface changes its topology of under the influ- 
ence of external pressure. Lisfhitz's original idea was 
extended to continuum theories of interacting fermions 
by Volovik 14] in the context of 3 He, where topologi- 
cal invariants were constructed for topological changes 
in the quasiparticle excitation spectrum, but thermody- 
namic signatures were not analyzed. Here, we generalized 
their pioneering work to the case interacting fermions in a 
lattice, and described thermodynamic and transport sig- 
natures of topological quantum phase transitions in the 
superfluid phase. 

Summary: We studied topological quantum phase 
transitions of ultracold fermions in optical lattices. For 
brevity, we focused on quasi-one-dimensional optical lat- 
tices, where the nature of the topological transitions in 
the superfluid state is most dramatic. We classified the 
quantum superfluid phases in accordance to their Fermi 
surface topologies and quasiparticle excitation spectrum. 
We showed that for s-wave superfluids there is no phase 
transition, but for p-wave (d-wave or /-wave) superflu- 
ids, quantum phase transitions of topological nature can 
occur, in particular from a phase with gapless quasipar- 
ticle excitations (BCS regime) to a phase of with gap- 
ful quasiparticle excitations (BEC regime). Finally, we 
showed that non-analytic behavior of the momentum dis- 
tribution, compressibility, sound velocities and superfluid 
density tensor as a function of interaction strength char- 
acterize well the topological phase transitions for p x , p y 
and p z , and f xyz pairing symmetries. We would like to 
thank NSF for support (Grant No. DMR-0709584) 
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